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It is shown that two gravitating scalar fields may form a thick brane in 5D spacetime. The 
necessary condition for the existence of such a regular solution is that the scalar fields potential 
must have local and global minima. 
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I. INTRODUCTION 

i> : 

' In recent years there has been a revived interest in theories having a greater number of spatial dimensions than the 
, three that are observed. In contrast to the original Kaluza-Klein theories of extra dimensions, the recent incarnations 
of extra dimensional theories allow the extra dimensions to be large and even infinite in size (in the original Kaluza- 
Klein theories the extra dimensions were curled up or compactified to the experimentally unobservable small size 
of the Planck length: 10^'^'^ cm). These new extra dimensional theories have opened up new avenues to explaining 
some of the open questions in particle physics (the hierarchy problem, nature of the electro-weak symmetry breaking, 
explanation of the family structure) and astrophysics (the nature of dark matter, the nature of dark energy) [l[ - 
In addition they predict new experimentally measurable phenomenon in high precision gravity experiments, particle 
. accelerators, and in astronomical observations. 
J> ' Most of the brane world models use infinitely thin branes with delta-like localization of matter. However these 
' models are generally regarded as an approximation since any fundamental underlying theory, such as quantum gravity 
or string theory, must contain a fundamental length beyond which a classical space-time description is impossible. It 
is therefore necessary to justify the infinitely thin brane approximation as a well-defined limit of a smooth structure 
- a thick brane - obtainable as a solution to coupled gravitational and matter field equations. One early example of 
a thick brane comes from the 5 dimensional model considered in [3| where one had a topologically non-trivial field 
. configuration for the scalar field. 
' In Ref. 0, the picture is presented that our universe is a dynamically localized 3-brane in a higher dimensional 
space ("brane world" ). As an example, the dynamics of the Nielsen- Olesen vortex type in six dimensional spacetime 
I ■ is adopted to localize our space-time within a 3-brane. At low energies, everything is trapped in the 3-brane, and the 
^jTj Einstein gravity is induced through the fiuctuations of the 3-brane. 

• • It is therefore of great interest to formulate some general requirements on the brane world design leading to the 

. !^ appearance of stable, thick branes having a well-defined zero thickness limit and able to trap ordinary matter. Taking 
, a physically reasonable stress-energy tensor it was shown that in 6 dimensions [6] and also higher extra dimensions 
^ ' one can trap all the Standard Model fields using gravity alone. 
. 5^ , In Ref's 0] thick brane world models are studied as Z2-symmetric domain walls supported by a scalar field with 
an arbitrary potential V{4>) in 5D general relativity and it was shown that in the framework of 5D gravity, a globally 
regular thick brane always has an anti-de Sitter asymptotic and is only possible if the scalar field potential V{(j)) has 
an alternating sign 

In Ref's 0] - 14 1 some properties of brane models was investigated: localization of gravity, graviton ground state. 



o 
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stability and so on. 

In Ref. a comparative analysis of localization of 4D gravity on a non Z2-symmetric scalar thick brane in both 
5-dimensional Riemannian space time and pure geometric Weyl integrable manifold is presented. 

Multidimensional space-times with large extra dimensions turned out to be very useful when addressing several 
problems of the recent non-supersymmetric string model realization of the Standard Model at low energy with no 
extra massless matter fields . 

In Ref. it is shown that two interacting non-gravitating scalar fields with a non-trivial potential may have a 
regular spherically symmetric solution. This solution shows that one can avoid the Derrick's theorem 18] forbidding 
the existence of regular static solution in the spacetime with the dimension greater 2 for scalar fields if the potential 
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has a local minimum besides global one. This result allows us to assume that the inclusion of gravitation may not 
destroy the regularity of similar solutions in 5D spacetime. In Ref. [l9| one example of spherically symmetric solution 
with a gravitating scalar field is given but in contrast with the solution that will be presented here the potential of 
the scalar field in Ref. is negative. 

The goal of this investigation is to show that there exists a new kind of thick brane solutions that is different with 
thick brane solutions found in Ref's [l^l [lH. We will show that the asymptotical behavior of one scalar field allow 
us to offer trapping of Maxwell electrodynamics and spinor fields on the brane. Especially it is necessary to note that 
the consideration of two scalar fields allow us to obtain the regular thick brane solution with the potential bounded 
from below. 



II. INITIAL EQUATIONS 



We consider 5D gravity + two interacting fields. The key for the existence of a regular solution here is that the 
scalar fields potential have to have local and global minima, and at the infinity the scalar fields tend to a local but 
not to global minimum. 

The 5D metric is 



ds^ — a{y)rif^i,dx'^dx" — dy^, 



(1) 



where — 0, 1, 2, 3; y is the 5*'' coordinate; jy^^ = { + 1, —1, —1, —1} is the 4D Minkowski metric. The Lagrangian 
for scalar fields (j) and x is 



(2) 



where A = 0, 1, 2, 3, 5. The potential V{4), x) is 



(3) 



where Vq is a constant which can be considered as a 5D cosmological constant A. We consider the case when the 
functions 0,% are (/)(?/), x(i/). The 5D Einstein and scalar field equations are 



dV{cl>,x) 
d(t> '■ 

dV{<l>,x) 
dx 



(4) 
(5) 
(6) 



where m is the 5D gravitational constant; Gab is the 5D metric ([T]) and G is the corresponding determinant. After 
substituting metric ([T]) into Eq's ([4|) - ([6]) we have the following equations 



x" + 4-x' 



1 

1 



12 , -^1 



,12 I 12 , '^l I i2 2\2 

b + X + Y (0 - m 



A2 



'2 / 2 



+ y (x^ - -2)' + 2^V- 2^0 



[2x^ + Ai ^^^/^ — , 
X [202 + A2 (x' - ml)] 



(7) 

(8) 

(9) 
(10) 



■ )'. Let us introduce the following dimensionless functions a/ y/>c/6 — > a, ipyGiJS — > 0, X\/ X; 



where ^^t-^ = 

ay 

2{m:/6) Vq — > Vq, mi_2\/ >f/3 mi^2, Ai,2/2 — > Ai,2 and the dimensionless variable y/\/x/& y. 
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After algebraical transformations Eq's ([7]) - (fTO]) have the following form 



6'' + x'') 



112 I a 



,2 2\2 A2 



a 

4-x' 
a 



+ A2 (x' 



(11) 

(12) 

(13) 
(14) 



It is easy to see that Eq. (|lip is the consequence of Eq. (IT^ : if we take a derivative from the LHS and RHS of Eq. 
then we shall receive Eq. (|lip . The boundary conditions are 



a(0) = ao, 

a'(0) = 0, 

0(0) = 00, 0'(O) = 0, 

X(0) = Xo, x'(0) - 0. 

The boundary condition pBjl - P^ and Eq. give us the following constraint 



(15) 
(16) 
(17) 
(18) 



Ai 



- jJ2 2 



(19) 



III. NUMERICAL INVESTIGATION 



For the numerical calculations we choose the following parameters values 

ao = 00 - 1, Xo - %/(I6, Ai = 0.1, A2 - 1.0. (20) 

We apply the methods of step by step approximation for finding of numerical solutions using the MATHEMATICA 
package (the details of similar calculations can be found in Ref. IITI , the corresponding MATHEMATICA program 
can be found in *.tar.gz file of the archived version of this paper 22]). 

Step 1. On the first step we solve Eq. p4[) (having zero approximations ao(y) = ao,Xo(2/) ~ mitanhy). The 
regular solution exists for a special value rri{ ^ only. For mi < m| ^ the function Xiiv) ~* +00 ^-nd for mi > m\ ^ 
the function Xiiu) —00 (here the index i is the approximation number). One can say that in this case we solve a 
non-linear eigenvalue problem: Xiiv) is the eigenstate and ^ is the eigenvalue on this Step. 

Step 2. On the second step we solve Eq. (|13p using zero approximation ao(y) for the function a{y) and the 
first approximation Xiiv) for the function xiu) from the Step 1. For 7712 < rrij i the function 4>i{y) — > +00 and for 
m2 > m2 1 the function 0i(j/) — > — cxd. Again we have a non-linear eigenvalue problem for the function 0i(y) and 

^24- 

Step 3. On the third step we repeat the first two steps that to have the good convergent sequence (j)* {y) , x* (y) ■ 
Practically we have made three approximations. 

Step 4. On the next step we solve Eq. (|lip which gives us the function ai(y). 

Step 5. On this step we repeat Steps 1-4 necessary number of times that to have the necessary accuracy of definition 
of the functions a* (y) , c/)* (y) , x* (y) ■ 

After Step 5 we have the solution presented on Fig. [1] These numerical calculations give us the eigenvalues 
m* « 2.122645756, « 1.3721439906788 and eigenstates a* (y) , cj)* (y) , x* (y) ■ The derived solution was verified 
by using the standard numerical method of solving the differential equations in the MATHEMATICA package (the 
corresponding MATHEMATICA program can be found in *.tar.gz file of the archived version of this paper [231). 

It easy to see that the asymptotical behavior of the solution is 

a{y) « a^e-'^^y, kl = \ [v, - ^m^^ , (21) 

dp{y) « mi + 0ooe-'=*^ = 2fc, + ^Akl + 2\iml, (22) 
X(y) « Xooe-'=-^ k^ = 2ka + ^Akl + ml-X2ml, (23) 



FIG. 1: The functions a* (y) , (p* (y) , x* (y) 



FIG. 2: The profile of dimensionless energy density. 



where aooi(t>ac ,X<x> are constants. The dimensionless energy density is 



1 



- 2Vo 



(24) 



and it is presented in Fig. [2l 

Taking into account that the quantity V ((/)(oo), x(oo)) is absolutely similar to a 5D cosmological constant, we can 
introduce a dimensionless brane tension 



a = 2 



3(2;)-y(0M,x(«3)) 



dy « 0.74. 



(25) 



According to Eq. (|22|) ([23]) one can define the thickness A of the presented thick brane as 

A w max {k^, k^} . 



(26) 



The key role for understanding why such regular solution may exist belongs to the fact that the potential ^ has the 
local and global minima. The profile of the potential V{(j),x) is presented in Fig. [3l 



IV. TRAPPING OF THE MATTER 



Now we would like to consider trapping of the electromagnetic and spinor fields on the above derived thick brane. 
The Lagrangian of interacting electromagnetic and scalar fields is taken from : 

Leff = -^FbcF^^ + u^'^AbA^ - n?ABA^, (27) 

where Fbc = Ab,c ~ Ac^b is the 5D electromagnetic tensor with 5-dimensional vector potential Ab{x^) and scalar 
field (p{y) depending only on the extra coordinate y; a - an arbitrary constant, m is the mass of vector field Ab- 
The 5D Maxwell equations will be: 

DcF^^ ^A^ix'')[a^' -m^]. (28) 

Let us rewrite Eq. (|28p as follows: 

D.F^'' + D5F^'> ^ A^{x^)[a(f>'' -m^]. (29) 
We will use the gauge A^ — and search for a solution of (|29p in the form: 

D^F'"' = 0, (30) 



chi 



-2 

phi 



FIG. 3: The profile of tfie potential 1/(0, x)- 



For the solution wc will use the following ansatz 

i^(a;^)=A^(x^)/(y), (32) 

where A-^{x'^) is the 4D electromagnatic potential function only on 4D coordinates. Then from Eq's ([30]) (|3T|) we will 
have 

D^F'"' = 0, (33) 

where the first equation is the usual 4D Maxwell equations on the brane. The solution of the second equation on 
the background of the thick brane is presented in FigUl Here it is necessary to note that again the regular solution 
f{y) exists for an exceptional value of the parameter m only. It is easy to see from Eq. (j34p : this equation is exactly 
Schrodinger equation with the potential acj)^ (which is a hole). Eq. (|34p has a regular solution describing a particle 
in a hole for an exceptional value of m that is an eigenvalue of the Schrodinger equation p4p . 

As one can see, the EM field is trapped on the 4D brane. In this case the electromagnetic fields in the bulk are 

Aix^) = Aix^fiy) (35) 

where f{y) is the exponentially decreasing function. 

Let us consider further the question about trapping of fermion fields on the brane. In the simplest case such a 
possibility was pointed out in Ref. psj at consideration of the brane model as the model of domain wall. In this work 
the model of one real scalar field 4> with two degenerated minima was introduced for description of the domain wall 
in 5D spacetime. In this case existing kink solution has its asymptotes in these minima with constant values of the 
field (f>. In our case similar situation occurs: two scalar fields (j), x create the system with two local minima, and the 
solutions tend asymptotically to one of these minima where the field x tends to zero and (p to the constant values as 
in the case from Ref. (23| . 

It allows us to investigate trapping of fermions on the brane for our case by analogy with Ref. [23 |. The curved 
space 5D gamma matrices are 



(36) 



FIG. 4: The function f{r) with a = 2 and m 1.6256513943551 



where 7^ and 7^ are the usual Dirac matrices in 4D theory 



/ 
/ 



a 

-3 



-1 
1 



(37) 



where cr' are usual Pauli matrices in flat spacetime. Then using the action for interacting scalar ip and fermion ^ 
fields we have 

= J d'^xdr{i^r^DA'i' ~ h0^) (38) 

here is the scalar field from the Lagrangian ([2|) and h is a constant. The Dirac equation can be written in the form 

iT^D^* - /i<?!)(y)* = 0. (39) 
Here — dA + T^, where pseudo-connection can be defined as follows [l^l 



where the vielbein is defined via gAB — ^a^b^mnj ^^'^ the inverse vielbein e^^ via g^^ = e^^e^r;*^^. 
case Ta = (0, 0, 0, 0, a'/a). 
Let us consider ansatz 

*(.T^) = V(^^)*o(y) 



For our 



(40) 



If we are interested in localization of zero modes, then, as it was shown in Ref. |25j . there are the solutions of Eq. 
(IM)) with 4D mass m = 0. For the zero mode j^Dfj^tp = 0, and the Dirac equation (1551) turns out in the equation: 



(41) 



where ' means the derivative with respect to r. Eq. (|¥T|) with account of P0)l has the following solution: 



5* — exp 



dr' i— + h(j){r'' 
a 



V(a;^), 



(42) 



where ip{x^) is the usual solution of 4D Weyl equation, and the condition 75 4*0 = —4*0 is taken into account. As 
it was shown in Section [1111 the sum (a'/a + hcf) tends asymptotically to some constant. So the zero mode (|^^ is 
localized near r = 0, i.e. on the brane, and decreases exponentially at large r: oc exp (— ms \y\). 
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Let us note that one can include the function x in Eq's ((27)) and ([38)) by the following way: 
a(t>'^ABA^ change ^ ^^^B and '^^^^ /i (0 + -^^ "^v]/ but it does not matter because the asymptotical 

behavior of the function is important only (as x '-^^ 0). 

Let us note that we do not consider trapping of scalar fields on this brane. The reason is very simple: we have 
shown exactly that two scalar fields with Lagrangian ([2]) are confined on the brane. The situation is even better: 
these scalar fields create the brane ! It is necessary note that in the process of numerical calculation we have obtained 
a domain wall solution without gravity, i.e. two scalar fields can create the solution with the planar symmetry and 
switching on the gravity does not destroy this solution. 



V. DISCUSSION AND CONCLUSIONS 

Now we would like to list the essential specialities of the presented solution: 

1. The existence of the solution crucially depends on the number of interacting scalar fields (n > 1) and the presence 
of the non-trivial potential V{(j), x) which has local and global minima. At the infinity the scalar fields tend to 
local minimum and the potential has alternating sign when r g [0, oo\ that to the existence of the presented 
solution. The numerical investigation shows that in the presence of one scalar field the similar solution does not 
exist. 

2. The advantage of the presented solution is that the asymptotical behavior (|^^ of the scalar field (j) allow us to 
obtain trapping of electromagnetic and spinor fields on the brane. 

3. Let us note that the thick brane solution presented here differs from the thick brane solutions presented in Ref 's 
and ^ that: 

(a) thick brane solution from Ref. [20| is obtained for the scalar field with the potential unbounded from below 
that in contrast with our potential ([3]) which is bounded from below. 

(b) in Ref. [2l[ the thick brane solution is obtained for scalar fields having non-trivial asymptotical topological 
structure in contrast with our solution. 

4. The solution is topologically trivial. It means that at the infinity two scalar fields do not form a hedgehog 
configuration in contrast with the thick brane solutions presented in Ref. (2l| . 

5. The quantity T^(0(oo), x(oo)) can be considered as a 5D cosmological constant A. 

6. In Ref. pJ7| it is shown that after some simplification and assumtions the SU(3) gauge Lagrangian can be 
reduced to the Lagrangian (U) describing interacting scalar fields cj) and x- This remark allows us to assume 
that a real thick brane can be formed by a 5D gauge condensate which is described by interacting scalar fields. 

7. According to the previous item ([6]) the 5D mechanism of trapping the matter on a thick brane may be similar 
to the confinement mechanism in 4D quantum chromodynamics. In this case trapping of the corresponding 
quantum gauge fields on the thick brane is non-perturbative and can not be investigated using Feynman diagram 
technique. 

8. If the thick brane is formed with the help of a gauge condensate then the problem of the stability of the thick 
brane becomes very non-trivial. It occurs because a non-static condensate has to be described in much more 
complicated manner than static condensate in Ref. It is connected to that fact that the change in time of 
quantum object is connected not only to the change of this quantity but also to the change of a wave function 
as well. 

9. From the mathematical point of view the presented solution is an eigenstate for a nonlinear eigenvalue problem. 
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